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We consider a family of unstructured problems, for which we propose a method for constructing 
analog, continuous-time quantum algorithms that are more efficient than their classical counterparts. 
In this family of problems, which we refer to as 'scrambled output' problems, one has to find a 
minimum-cost configuration of a given integer-valued n-bit function whose output values have been 
scrambled in some arbitrary way. Special cases within this set of problems are Grover's search 
problem of finding a marked item in an unstructured database, certain random energy models, 
and the functions of the Deutsch-Josza problem. We consider a couple of examples in detail. In 
the first, we provide a deterministic analog quantum algorithm to solve the seminal problem of 
Deutsch and Josza, in which one has to determine whether an n-bit boolean function is constant 
(gives on all inputs or 1 on all inputs) or balanced (returns on half the input states and 1 on 
the other half). We also study one variant of the random energy model, and show that, as one 
might expect, its minimum energy configuration can be found quadratically faster with a quantum 
adiabatic algorithm than with classical algorithms. 
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I. INTRODUCTION 

The paradigm of Adiabatic Quantum Computation 
(AQC) proposed by Farhi et al. [l| about a decade ago 
is a simple yet intriguing approach to problem solving 
on a quantum computer. Unlike the leading paradigm 
of circuit-based quantum computing, AQC is an analog 
continuous-time method that does not require the design 
and use of quantum gates. As such, it can in many ways 
be thought of as a simpler and perhaps more profound 
method for performing quantum computations that is 
also easier to implement experimentally 0, Q . 

Even though AQC has been shovifn to bepolynomially- 
equivalent to circuit-based computation [j, |5|, and de- 
spite intensive research in the area (see, e.g., Refs. 
Il4| and references therein), to date, there are almost 
no clear-cut concrete examples for efficient quantum- 
adiabatic algorithms that reveal the potentially-povifcrful 
"fully-quantum" capabilities encompassed in AQC. One 
possible reason for that is, presumably, that there is usu- 
ally no obvious way to 'tailor' the adiabatic algorithm to 
the specific problem being examined, and to make use of 
the structure of the problem to, for example, modify the 
beginning Hamiltonian in a clever way that would speed 
up the computation (a notable exception is Ref. jl5l|). 
For most of the interesting optimization problems, be- 
ing able to do so, may be as hard as solving the original 
problem itself [301 ■ 

Interestingly, one of the few problems for which quan- 
tum speedup has been obtained in the context of AQC, is 
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Grover's unstructured search problem [16|, in which one 
searches for a marked item in an unstructured database. 
Roland and Cerf [T^] have demonstrated that while the 
application of the adiabatic algorithm to Grover's prob- 
lem with a linear rate results in a running time that is of 
order N - N being the number of items in the database 
- a carefully chosen variable rate of the adiabatic pa- 
rameter yields a running time that scales like \/N, i.e., 
a quadratic speed-up is gained, similarly to the original 
result by Grover found for the circuit-based model |3l| . 

The adiabatic algorithm for Grover's problem utilizes 
the concept of local adiabatic evolution, in the frame- 
work of which the adiabatic parameter is varied not at a 
constant rate but rather at a variable rate, slowing down 
in the vicinity of the minimum gap and speeding-up in 
places where the gap is large. 

Local adiabatic evolution however can only be effi- 
ciently used in cases where one has proper knowledge 
of the exact behavior of the gap and relevant matrix ele- 
ments of the system for the problem in question. This is 
normally not the case. In the Grover problem, the ability 
to compute the gap and matrix element of the problem 
stems from prior knowledge of the spectrum of the prob- 
lem Hamiltonian, which, ultimately, reduces the problem 
into a simple two- level system (l7| . 

In what follows, we consider a family of unstructured 
problems, which we refer to as 'scrambled output' mod- 
els, and show how one may utilize knowledge of the spec- 
trum of the Hamiltonian of the problem to find analog, 
continuous-time, algorithms that are more efficient than 
their classical analogues. In that sense, this family of 
problems is a generalization, or an extension, to the prob- 
lem solved by Roland and Cerf. 

In scrambled output problems, one has to find a mini- 
mum input configuration (i.e., argmin) of an n-bit func- 
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tion whose set of outputs (and their multiphcities) is 
given in advance, up to an unknown constant offset. The 
exact mapping between tlic = 2" input configurations 
and the various outputs is also not given (i.e., it is as 
though the outputs of a known function have been scram- 
bled in some arbitrary way). As we also discuss later, 
special cases in this family of problems are the unstruc- 
tured database search problem considered by Roland and 
Cerf ITj, certain variants of the random energy model 
and the functions of the Deutsch-Josza problem. 

We illustrate the manner in which AQC may be used 
to solve scrambled output problems by constructing effi- 
cient analog continuous-time algorithms for two specific 
examples: the Deutsch-Josza [l9| problem for which we 
find an efficient 0(f) deterministic solution, and a vari- 
ant of the random energy model, for which the minimum 
energy configuration is found quadratically faster than 
the corresponding classical algorithms. 

The paper is organized as follows. In the next section 
we briefly discuss the principles of the Quantum Adi- 
abatic Algorithm that is the heart of AQC, and with 
which the above models are solved. In Sec. Illli we de- 
scribe scrambled output problems in detail. We then 
study two examples. In Sec IIVI we suggest an analog al- 
gorithm for the Deutsch-Josza problem, and in Sec.lVlwe 
consider a variant of the random energy model. Finally, 
we conclude with a few comments in Sec. IVIl 



For adiabatic processes, the adiabatic profile function 
s{t) must be chosen such that the evolution of the system 
is slow. In the simple case where s(t) varies from zero to 
one at a constant rate, the runtime T must be chosen 
to be large enough so that the adiabatic approximation 
holds: this condition determines the efficiency, or com- 
plexity, of the QAA. A condition on T can be given in 
terms of the instantaneous eigenstates ||?7t)| and eigen- 
values {E„i} of the Hamiltonian H{s), as [3, 



^ _ 1 maxs Vbi(s) 
e mills g^{s) 



(2) 



where g{s) is the first excitation gap Ei{s) — Eq{s) and 
Voi{s) = |(0|diJ/ds|I)| (in our units h = 1). Here, e 
is a small number inversely proportional to the running 
time of the algorithm. The smaller e is chosen to be, the 
slower the evolution of the system will be and the larger 
the probability of success will become. For an adiabatic 
process to have zero error, the running time must tend 
to infinity. 

As discussed in the Introduction, one could gain sig- 
nificant speed-up in calculation times in cases where the 
concept of local adiabatic evolution can be utilized. 
This is done by formulating a 'local' Landau-Zener con- 
dition for each value of the adiabatic profile function s(t), 
namely: 



II. QUANTUM ADIABATIC ALGORITHM 
(QAA) 

The Quantum Adiabatic Algorithm (QAA) provides 
the general approach for solving optimization problems 
on an analog continuous-time quantum computer [ij. 
Within the framework of the QAA, the solution to an 
optimization problem is encoded in the ground state of 
a Hamiltonian Hp. To find the solution, the QAA pre- 
scribes the following course of action. As a first step, the 
system is prepared in the ground state of another 'driver' 
Hamiltonian Hfi. The driver Hamiltonian is chosen such 
that it does not commute with the problem Hamiltonian 
and has a ground state that is fairly easy to prepare. 
As a next step, the Hamiltonian of the system is slowly 
modified from Hd to Hp, using the linear interpolation, 
i.e.. 



H{s) = sHp + {l-s)Hd, 



(1) 



where s{t) is a parameter varying smoothly with time, 
from at i = to I at the end of the algorithm, t = T- 
If this process is done slowly enough, the adiabatic theo- 
rem of Quantum Mechanics (see, e.g., Refs. [l^l and plj ) 
ensures that the system will stay close to the ground 
state of the instantaneous Hamiltonian throughout the 
evolution, so that one finally obtains a state close to the 
ground state of Hp. At this point, measuring the state 
will give the solution of the original problem with high 
probability. 



ds 
dt 



< e 



(3) 



where e. 17(5) and Vbi(s) are as in Eq. ([25]) . 

Before moving on to describe scrambled output prob- 
lems, we note here that in special cases, the 'machinery' 
of the QAA may also be used differently than described 
above, namely, in a not-necessarily adiabatic fashion. By 
'non-adiabatic' use, it is meant that the state of the sys- 
tem may 'wander away', at least to some extent, from the 
instantaneous ground-state of the system in the course of 
the algorithm before returning to it towards the end of 
the run. This concept is discussed later in greater detail. 
One known example for such 'non-adiabatic' computa- 
tion is the analog version of Grover's search algorithm 
suggested by Farhi and Gutmann in Ref. [Ill ■ 



III. SCRAMBLED OUTPUT PROBLEMS 

In what follows, we consider a family of problems, that 
we refer to as 'scrambled output' problems, and which, 
as we illustrate next, can be solved more efficiently on an 
analog quantum computer than on a classical computer. 
As we explain in what follows, this family of problems 
is a generalization, or an extension, of the unstructured 
database search considered by Roland and Cerf |I7|. 

In 'scrambled output' problems, one is asked to find a 
minimizing input configuration of an integer- valued n-bit 
function whose A' + I output values fo < fi < ■ ■ . < /k 
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and their multiplicities mo, mi, . . . , mx are known in ad- 
vance, up to an unknown constant offset, denoted by eg, 
that is added to all outputs. Here, we shall assume that 
the number of distinct output values, if + 1 is much 
smaller than the 2" actual eigenvalues. While the out- 
put values of the function arc known, it is not known 
which output belongs to which input bit-configuration; 
the output values arc 'scrambled' in some arbitrary way. 
The additional unknown constant offset serves to further 
complicate the problem: It makes the minimum-energy 
configuration harder to find for classical algorithms, but, 
as we shall see, will have no effect on the quantum ones. 

Special cases in this family of problems are: i) The 
original Grover's search problem considered by Roland 
and Cerf [TtI in which /o = with multiplicity mo = M 
(where M is the number of marked items or 'targets') and 
/i = 1 with mi = N — M . ii) The Deutsch-Josza input- 
functions (which we discuss in the next section) for which 
both the balanced and constant functions are described 
by /o = and /i = 1, with multiplicities mo = mi = N/2 
(balanced) and mp/i = and mi/Q = N (constant), iii) 
A variant of the 'random energy model' for which fj ~ j 
with j ~ 0..n and mj = (") (n-choose-j). This model 
has been previously considered in Ref. [2^ in the context 
of the QAA. 

Since the output values of the given function are scram- 
bled in an arbitrary way, the typical running time of any 
classical search algorithm designed to find a minimum- 
cost configuration will scale as 0{N) (unless of course 
there are exponentially many minimizing configurations) . 
For example, a classical deterministic solution to the 
problem will require in the worst case N — mo evalua- 
tions of the scrambled-output function. 

In quantum-algorithmic terminology, the problem de- 
scribed above can be encoded in a matrix F (which we 
take to be diagonal in the computational basis), whose 
eigenvalues are known but have been shuffled around and 
added an unknown offset bq. We will take this F to be 
the problem Hamiltonian of a QAA procedure: 

Hp = eo (4) 
+ diag[7r( /o, . ^. , /q , /i, . ^. , /i , . . . , /k, . , Jk) ] , 

mo times rai times rriK times 

where 7r(-) is an arbitrary permutation of its list of argu- 
ments. 

The lack of structure exhibited by the problem Hamil- 
tonian above does not seem to allow for any clever 'tailor- 
ing' of QAA to the problem, as no information about the 
problem can be incorporated into the driver Hamiltonian 
that would help in speeding up the computation. There- 
fore, a driver Hamiltonian that treats all computational- 
basis eigenvectors in exactly the same way will presum- 
ably be an optimal choice. We thus take as driver Hamil- 
tonian, the one-dimensional projection onto the equal- 
superposition state: 

Hd^^Eo\(f>){^\, (5) 



where Eq is a positive constant that provides a scale to 
the driver Hamiltonian and \(f>) is the fully-symmetric 
state 

1 ^ 

The ground-state of this driver Hamiltonian, i.e., the 
state that the system is prepared in, for the adiabatic 
evolution, is 

m = 0)) = , (7) 

whose energy is ^Eq. 

We note here that both the scrambled problem Hamil- 
tonian and the above driver Hamiltonian have been 
shown to pose severe limitations on the efficiency of any 
QAA constructed using either of the two [l^ . We shall 
address this matter later on. Here, however, the above 
choice of driver Hamiltonian has an attractive prop- 
erty, at least as far as scrambled problem Hamiltonians 
are concerned. Its symmetry makes it invariant under 
any permutation of the eigenstates of the computational 
(problem Hamiltonian) basis. Specifically, it is invariant 
under the (unknown) unitary transformation that 'un- 
scrambles', i.e., orders, the eigenvalues of the problem 
Hamiltonian. It thus follows that the Schrodinger equa- 
tion of the QAA for this type of problems can be writ- 
ten down explicitly, and subsequently solved, at least in 
principle. This property of the system can then be used, 
as we show next, to find a (not necessarily) adiabatic 
path s{t) such that at the end of the algorithm, the final 
state of the system lies very close to (or in some cases, is 
precisely) the ground-state of the problem Hamiltonian. 
Once an optimal path is found, a measurement of the 
energy at the end of the run, will immediately reveal the 
minimum energy eo -l- /o (and therefore also the value 
of the unknown offset eo), along with the corresponding 
minimum-energy configuration. 

Explicitly, the Schrodinger equation for a QAA on a 
scrambled output problem, i-^\ip{t)) = H\ip{t)), simpli- 
fies to the (K + l) coupled first-order ordinary differential 
equations: 

K 

icj ^ -{1 - s)^rikCk + sfjCj , (8) 

where j — 0..n. Here, rjj = mj/N are the ratios of the 
multiplicities of occurrences of the j-th eigenvalue fj to 
the total number of eigenvalues N. The various Cj are 
the {K + 1) distinct wave-function amplitudes that cor- 
respond to the different values fj. We have thus ended 
up with a compact set of equations which needs to be 
solved or analyzed in some meaningful way. 

A few remarks are now in order: i) Note that the ef- 
fective dimension of the problem is K + l, i.e., precisely 
the number of distinct /j's. ii) In terms of the differ- 
ent amplitudes Cj, the initial-state of the system corre- 
sponds to Cj{t = 0) = 1 for all j, where the normal- 
ization chosen here is X^jLo % I'^j P ~ 1- ^otc that the 
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evolving wave- function depends on iV, the size of the sys- 
tem, only through the various ratios, iii) The addition 
of an unknown offset eq to the various fj only affects the 
(immaterial) global phase of the wave-function, and may 
therefore be completely removed from the equations. 

In what follows, we utilize the ideas presented above, 
to obtain efficient analog algorithms for two exem- 
plary problems for which no efficient classical solution is 
known. We first consider the famous Deutsch- Josza prob- 
lem [l9l | , for which we find an efficient constant- runtime 
deterministic solution, that stands in contrast to the ex- 
ponential complexity of the corresponding deterministic 
classical algorithm, and is on par with the circuit-based 
quantum Deutsch- Josza algorithm. We next consider a 
variant of the random energy model already mentioned in 
the previous section, for which the minimum energy con- 
figuration is found quadratically faster than on a classical 
computer. 



IV. AN EFFICIENT DETERMINISTIC 
ALGORITHM FOR THE DEUTSCH-JOSZA 
PROBLEM 

In the Deutsch- Josza problem, we are given a black 
box quantum computer (an oracle) that implements a 
boolean function of n bits that is either constant, i.e., 
gives on all inputs or 1 on all inputs, or balanced, i.e., 
returns 1 for exactly half of the input states and for the 
other half. The task then is to determine with absolute 
certainty whether the function is constant or balanced 
with as few calls as possible to the oracle. 

A classical deterministic algorithm would require 
N/2 + 1 evaluations of the function in the worst case (and 
two evaluations in the best case). The seminal result by 
Deutsch and Josza [l^ was the construction of a deter- 
ministic circuit-based quantum algorithm that requires 
only two calls to the oracle (later reduced to only one 
call by Cleve et at. [l^), thus providing the first exam- 
ple of a quantum computation that is exponentially more 
efficient than the best corresponding classical algorithm. 

In what follows, we show that one can use the princi- 
ples discussed in the previous section to find an analog de- 
terministic quantum computation to solve the Deutsch- 
Josza problem with a runtime that is on par with the 
analogous circuit-based algorithm, i.e., a runtime that 
does not scale with input size, and so is exponentially 
faster than the corresponding deterministic classical al- 
gorithm. 

The algorithm we propose here is very similar to adia- 
batic algorithms of the usual AQC paradigm but with 
one important distinction. While with adiabatic pro- 
cesses the success of the algorithm depends on the state 
of the system being close at all times to the instantaneous 
ground state of the evolving Hamiltonian, here shall 're- 
lax' this condition. We will show that by allowing the 
state of the system to 'detach' itself at least to some ex- 
tent from the instantaneous ground state, one can in fact 



obtain a deterministic (i.e., zero probability of failure) 
non-adiabatic efficient algorithms, for the Deutsch- Josza 
and certain other problems (see, e.g., Ref. psj). 

A. The adiabatic Deutsch-Josza equations 

In the context of 'scrambled output' problems, the 
Deutsch- Josza oracle corresponds in the 'constant' case 
to having /o = or /o = 1 (with toq = 2" = iV), and 
/o = and /i = 1 with multiplicity toq — 'tii — N/2 in 
the balanced case. The task here is to be able to distin- 
guish between the two cases. In this problem, we shall 
assume, for simplicity, that the constant offset cq is zero 
(even though in principle it could take on any value). 

To find a deterministic efficient solution to the 
Deutsch-Josza problem, we construct a continuous-time 
algorithm that is optimized for the case where the input 
matrix F is balanced, for reasons that will become clear 
shortly. In the balanced case, the problem Hamiltonian 
has two distinct values and the wave function \ip{t)) has 
only two distinct components. It may therefore be writ- 
ten as 

m)) = coimo) + ciit)\A) , (9) 

where j-^o) and \ipi) are the equal superpositions of all so- 
lution (zero-energy) states and non-solution (energy one) 
states, respectively: 

where the N/2 states m G Ai are eigenstates of the com- 
putational basis with zero eigenvalues (all other states 
have eigenvalue one) . The Schrodinger equation, Eq. ^ , 
thus becomes the two coupled equations (here, we fix the 
energy scale of the driver Hamiltonian at Eq = 1): 

tco =-i(l-s)(co + ci) (12) 
ici = -i(l - s)(co + ci) +sci , (13) 

As previously mentioned, note that N, the size of the 
problem, does not appear in the equation. 

The time-dependence of cp and ci is easily transformed 
into an s dependence. The time-derivative transforms in 
the usual way: d/dt = (ds/di) x d/ds, and the two cou- 
pled equations now become one, with the new complex- 
valued dependent variable r(s) = ci(s)/co(s): 

2ir'{s) = t'{s) [-1 + S + 2sr{s) + (1 - s)r{s)^] , (14) 

with the initial condition of r(0) = 1 corresponding to 
the fully-symmetric initial state. Here, the prime symbol 
(') stands for differentiation with respect to the new in- 
dependent variable s and t'{s) = 1/ (ds/di) encodes the 
non-adiabatic path. 
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As a next step, we split the above equation into its real 
and imaginary parts, re-expressing the complex-valued 
r(s) by real- valued components of magnitude and phase: 



V p(«) 



(15) 



Here, p{s) is the probability of the system to be in the 
symmetric zero-energy state {i/jq) at any given s, and 4>{s) 
is the phase, which is easily solved for, resulting in a 
single real- valued equation for the probability p{s). 

A deterministic algorithm whose end-state is precisely 
the equal superposition of the zcro-cncrgy states, IV'q), 
may be obtained if we satisfy the requirement that the fi- 
nal state of the system at s = 1 will be precisely p{l) ~ 1 
(equivalently, r(l) = 0). To make it so, we shall treat 
the Schrodinger equation, in this example Eq. (fT4)) . as 
an equation on the 'path' t'{s) rather than on the wave 
function itself [in our case, p(s)]. Put differently, we first 
choose the wave function (r(s) here) or rather, the prob- 
ability profile p(s), such that the probability of obtaining 
the desired solution is exactly one at the end of the run, 
and only then obtain the non-adiabatic path t'{s) con- 
sistently with the Schrodinger equation of the system, 

Eq. (m. 

For the Deutsch-Josza problem, an explicit expression 
for the path t'{s) as a function of the probability profile 
p{s) may be obtained: 



p'{s) 



Vp(s)(i-p(s))yr^^ 



where 



lis) 



i-.(i + 2p(.) + 2(i-.)/;^d.- 

2{l-s)y/p{s){l-p{s)) 



(16) 



(17) 



Presumably, there are infinitely many choices for p(s) 
that yield a desired path. One however must make cer- 
tain that p{s) satisfies the boundary conditions, p{0) = 
1/2 (corresponding to r(0) = 1), and p{l) — 1 (corre- 
sponding to r(l) = 0), and that the profile t'(s) is well- 
defined everywhere and yields a finite running time for 
the algorithm, namely, that T = /q i'(s)d.s < oo. One 
example for such a probability profile is 



P,(s) = i(l + 6.s2-8s3 



3^4) . 



which immediately leads to 



6^2 



Vs(l-s)[4-^9s(l-s)] 



(18) 



(19) 



Figure [T] shows the "adiabatic" parameter s* (i) for the 
above example, obtained by integration and inversion of 
t'^(s) above. The inset shows p*(t), the generating prob- 
ability profile. 
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FIG. 1: The path St(f) generated by the probability profile 
p*(s) given in Eq. (|18p . The probability p*(s) is plotted in 
the inset. 
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FIG. 2: Probability of being in the non-solution superposi- 
tion, 1 — p,(s), throughout the non-adiabatic evolution for the 
example discussed in the text (solid line) compared against 
the probability profile of the adiabatic case (dashed line). 
While the latter probability reaches at the end of the pro- 
cess a small yet strictly nonzero value at s = 1, in the former 
case the probability of failure drops precisely to zero. 



In Fig. [21 the probability of being in the non-solution 
superposition, namely, 1 — p*(s), is plotted as a function 
of s on a logarithmic scale. As can be seen in the figure, 
this probability (solid line) drops sharply to zero as s 
approaches one. This is in contrast with the probability 
of failure in the adiabatic case (dashed line) generated 
by the adiabatic path t' (s) = [e(l — 2s(l — s))]"^ which 
does not vanish as s = 1 (here, e = 0.001). 

As mentioned above, there arc presumably many paths 
s(t) that produce, for a carefully determined runtime, the 
desired end probability of one. We note here, that even 
the simple choice of s{t) = 1/2, while not conforming to 
the boundary conditions of s(0) = 1 and s(l) = 1, also 
produces the desired probability. The running time in 
this case, T = V^tt « 4.443, is however, substantially 
longer than the one found above. 
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B. The algorithm 

Now that it has been estabhshed that a deterministic 
algorithm to reach the zero-energy superposition for a 
balanced Deutsch-Josza problem Hamiltonian exists, we 
construct an efficient QAA for the Deutsch-Josza prob- 
lem as follows. Within the quantum-analog Deutsch- 
Josza algorithm, one executes the QAA twice: Once with 
a problem Hamiltonian of Hp = F and a second time 
with Hp = \ ~ F. Note that if F is balanced (constant) 
then (1 — F) is also balanced (constant). In both runs, 
the profile function to be used is one which guarantees 
the end-state to be the zero-energy superposition in the 
case where F is balanced (an example for such a path 
has been given above). 

Now, if F is constant, the problem Hamiltonian will 
have no effect on the state of the system in cither run, 
regardless of the profile function chosen for the run. The 
state will remain in the initial state throughout the 
evolution, except for an undetectable global phase. This 
is due to the constancy of the problem Hamiltonian. A 
measurement of the z-magnetization at the end of the 
run will thus pick out a random bit configuration that 
will have energy in one run and energy 1 in the other. 

If on the other hand F is balanced, the chosen path 
ensures us that at the end of each of the two runs, at 
s = 1, the final state of the system will be the equal 
superposition of the zero-energy states with probability 
1 . A measurement along the z-direction will thus produce 
in this case a zero-energy state in both runs. 

The resulting pair of energy measurements therefore 
distinguishes between a balanced function (both energy 
readouts are 0) and a constant function (energy read- 
outs of and 1), deterministically. Since t'{s) is iV- 
indcpendcnt. the running time is 0(1), i.e., it does not 
scale with the size of the problem. The above scheme thus 
provides a constant-runtime deterministic algorithm for 
solving the Deutsch-Josza problem. 

V. AN EFFICIENT ADIABATIC SOLUTION TO 
THE RANDOM ENERGY MODEL 

In the following example, we consider a different model, 
namely, a variant of the random energy model (26l . [27j , 
previously considered under somewhat different settings 
in Refs. [l^, [2^ in the context of the quantum adia- 
batic algorithm. Here, we show that one can utilize the 
knowledge of the spectrum of the problem Hamiltonian, 
to obtain an optimal adiabatic path that in turn yields 
a quantum adiabatic algorithm to find the minimum- 
energy configuration. The resulting algorithm is found 
to be, as expected, quadratically-faster than any classi- 
cal search algorithm. 

In this variant of the random energy model, the eigen- 
values of the problem Hamiltonian are integers in the 
range k G {0..n} where each value k appears precisely 
(^) times. The iV = 2" eigenvalues are distributed ran- 



domly along the diagonal. As discussed earlier, one could 
further complicate the problem by considering adding an 
arbitrary constant offset to each of the k values. Here 
too, one has to construct an efficient algorithm designed 
to find the minimum-energy bit configuration (i.e., the 
configuration that corresponds to A: — 0), and to de- 
termine how the running time of the algorithm scales 
with the dimension of the Hilbert space of the problem, 
N — 2". Note that for the corresponding classical algo- 
rithm, the average running time would scale like 0{N), 
as the ground state has multiplicity one and there is no 
efficient searching routine to find it, due to the lack of 
structure of the problem. 

For what follows, we choose the energy scale of the 
driver Hamiltonian to he Eq = n for computational con- 
venience. The Schrodinger equation, Eq. ([S]), becomes 
for this problem a set of only n + 1 linear first-order dif- 
ferential equations, corresponding to the n + 1 distinct 
energies of the problem Hamiltonian. The characteristic 
polynomial of the reduced Hamiltonian, as a function of 
the adiabatic parameter s, is simplified to the equation: 

fc=0 ^ ^ 

The solutions A to this equation are the eigenvalues of 
the Hamiltonian, and may be obtained analytically in 
the large n limit. 

Figure[3]shows the three lowest energy levels of the sys- 
tem as a function of s for n = 20 (the energy levels are 
qualitatively similar for larger all n values as well). The 
two avoided crossings between these three levels com- 
pletely determine the behavior of the gap and the matrix 
element of the system for all s values. A simple pertur- 
bation analysis of the system reveals that one needs only 
to consider the projection of the Hamiltonian into the 
subspace spanned by i) the ground state of the driver 
Hamiltonian, ii) the ground state of the problem Hamil- 
tonian and iii) the equal superposition of the first n ex- 
cited states of the problem Hamiltonian (all of which have 
the same energy). Analysis of the three- level system in 
the large n limit enables us to obtain rather easily the 
large-n asymptotic behavior of the gap: 

(n-^s <.<3^ 
g(s) = \'-^s-n 3^<,<_^ (21) 

In the immediate vicinity of the minimum gap, we find 
that 

ming(s) = -2-"/^ (22) 
s 6 

where the critical point is approaches =2/3 for large 
n. 

While the matrix element Vqi (s) may also be computed 
using perturbation theory, it is easier to bound it from 
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.... E2(S) X 
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0.0 0.2 0.4 0.6 0.8 1.0 



s 

FIG. 3: The three lowest energy levels for the n = 20 random 
energy model considered in the text. The inset is a blowup of 
the region of the of avoided crossings between the levels. 

above by simpler considerations: 




10 20 30 40 50 60 
n=log2N 

FIG. 4: The running time eT of the random energy model adi- 
abatic algorithm as a function of n, the number of input-bits, 
on a logarithmic scale. As the figure indicates, the running 
time scales like 



l/oi(s) = |(0|di//ds|l)| 

= -l(o|^.|i)l = 



(23) 



1 - s' 



l^^pll) 



The above matrix elements for both the problem and 
driver Hamiltonian arc easily bounded by their matrix 
norms to give: 



Vbi(s) < min 



1 - s 



= 2n. 



(24) 



While the above bound is by no means tight (the ac- 
tual tight bound is also linear in n albeit with a smaller 
constant), the above bound will suffice for our purpose. 

A constant-rate adiabatic scheme, following 
Eq. ([^5)) . will yield QAA running times that are 
0(maxs Voi(s)[emins 5(s)^]~^) = 0{N/ne), and would 
therefore provide no real advantage over the performance 
of classical algorithms. 

Nonetheless, our ability to calculate the gap of the sys- 
tem accurately, along with proper bounds on the matrix 
element, allows for the construction of a locally-adiabatic 
path that would in turn yield significantly shorter QAA 
running times. The bound on the matrix element, com- 
bined with the gap g{s) found in Eqs. (PT|) and , 
using the local adiabatic evolution condition, Eq. ([3|), 
allows us to construct the locally-adiabatic path which, 
upon integration, yields the following scaling on the run- 
ning time in the large n limit: 



r 



ds 



(25) 



- 0(maxVbi(s)[emin5(s)]-i) = 0(-^) . 
s s ne 

To verify that the large-scale analysis presented above is 
correct, we plot in Fig. 3] the actual running time as it 
was calculated by integrating the exact-numerical gaps 



and matrix elements for up to n = 60 (crosses, the e pre- 
factor has been removed) as a function of the number of 
input bits n, on a logarithmic scale. The solid line in 
the figure corresponds to eT = ^2"/^~^ = which 
turns out to be the true large-n behavior of the running 
time, and which provides a verification for the large n 
scaling analysis presented above. 

We have thus shown how to construct a quantum adia- 
batic algorithm with which the random energy model can 
be solved quadratically faster than any existing classical 
algorithm. 

We note here that the above analysis is in accord with 
theorems proved in Ref. [l^ by Farhi et al. which have 
shown that the algorithmic efficiency of the QAA with 
either a one-dimensional projection as the driver Hamil- 
tonian or an unstructured problem Hamiltonian (under 
certain conditions) is bounded by an 0{\/N) of the run- 
ning time. Here, we have used both. 



VI. SUMMARY AND CONCLUSIONS 

We have demonstrated how one can construct analog, 
continuous-time, quantum algorithms for a family of un- 
structured problems, namely, scrambled output. We have 
shown that for this type of problems, a compact set of 
dynamical equations may be written down explicitly, and 
subsequently solved or analyzed to enable finding opti- 
mized adiabatic or non-adiabatic paths, that yield algo- 
rithms that are more efficient than corresponding classi- 
cal algorithms. 

We considered two specific examples. In the first, a 
simple prescription for solving the Deutseh-Josza prob- 
lem was given. We have shown that the running time of 
the proposed analog algorithm does not scale with the 
size of the problem, i.e., that it is exponentially faster 
than the classical algorithm and on par with the circuit- 
based quantum algorithm result. 
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In a second example, we provided an adiabatic solution 
to the problem of finding the ground state of a random 
energy model, in which the eigenvalues are taken from 
a binomial distribution but are then shuffled around. In 
this example, the provided algorithm was found to be, 
as one might expect, quadratically faster than the corre- 
sponding classical one. 

It would be of interest to sec how this family of prob- 
lems may be further generalized, thereby expanding the 
scope of applicability of continuous-time quantum com- 
puting. We hope that the method presented here will 
open a way for other adiabatic as well as determinis- 
tic non-adiabatic efficient algorithms that would further 



demonstrate the power and potential encompassed in 
analog continuous-time quantum computation. Specifi- 
cally, we hope that the above algorithm will help in fur- 
ther pinpointing the precise equivalence between circuit- 
based and adiabatic quantum computing. 
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